NON-EXTINCTION OF A FLEMING- VIOT PARTICLE MODEL 
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Abstract. We consider a branching particle model in which particles move inside a Eu- 

0^ ' clidean domain according to the following rules. The particles move as independent Brown- 

O ■ 

■ ian motions until one of them hits the boundary. This particle is killed but another randomly 

chosen particle branches into two particles, to keep the population size constant. We prove 



in 
o 



that the particle population does not approach the boundary simultaneously in a finite 
time in some Lipschitz domains. This is used to prove a limit theorem for the empirical 
distribution of the particle family. 



1. Introduction 

^ . The paper is concerned with a branching particle system = {X^, . . . ,X^) in which 

individual particles move as independent Brownian motions and die when they hit 
the complement of a fixed domain D C M^. To keep the population size constant, whenever 
any particle X^ dies, another one is chosen uniformly from all particles inside D, and the 
■ chosen particle branches into two particles. Alternatively, the death/branching event can be 
^ ' viewed as a jump of the j-th particle. See Section [5] for a more detailed description of the 
construction. 

Let Tfc be the time of the k-th jump of X^. Since the distribution of the hitting time of 



^ ■ dD by Brownian motion has a continuous density, only one particle can hit dD at time r^, 
for every k, a.s. The construction of the process is elementary for all t < t^o = limfe^oo Tfc- 
^ i However, there is no obvious way to continue the process Xj after the time Tj, if Tqo < oo. 
. Hence, the question of the finiteness of Too is interesting. Theorem 1.1 in lOj asserts that 
Too = oo, a.s., for every domain D. Unfortunately, the proof of that theorem contains an 
irreparable error (see Example 15 . 31 below) . The cited theorem might be true but it appears to 
be much harder to prove that the original incorrect argument might have suggested. Example 
iven below shows that result cannot be generalized to arbitrary Markov processes. Lobus 
19|) recently proved that Tqo = oo, a.s., in Euclidean domains that satisfy the internal ball 
condition. Another argument showing that Too = oo, a.s., in domains satisfying the internal 



ball condition is implicit in the proof of Theorem 1.4 of lOj. 
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In this article, we will prove that t^o = oo, a.s., if the domain D C M'' is Lipschitz with 
a Lipschitz constant c depending on d and the number of particles — see Theorem 15.11 
and Remark 15.21 below. In addition, we prove theorems on existence and the form of the 
stationary distribution of the process X^, generalizing those in 10|] — see Section [71 

We use this attempt to rectify an error in an earlier paper to introduce two new techniques. 
In the end, these techniques may have greater interest or significance than the main theorems. 
The first technique, developed in Section HJ is the construction of a process of Brownian 
excursions in a cone, with all excursions starting at the vertex. Such a process exists only in 
cones with certain angles. The construction is combined with a coupling argument to provide 
a "lower bound" for Xj, in an appropriate sense. The process constructed from Brownian 
excursions is simpler to analyze than X^. 

The second technique is a new type of boundary Harnack principle (see Section E]). The 
standard boundary Harnack principle compares two functions satisfying a PDE with the 
same operator, for example, Laplacian, and different boundary conditions. Our new version 
of the boundary Harnack principle compares a harmonic function with a function u satisfying 
Au = — 1. The reason for proving the new form of the boundary Harnack principle is that 
it allows one to compare certain probabilities and expectations, and then use a method of 
proof that goes back at least to Davis [3]. The "new boundary Harnack principle" has 



been proved independently by Atar, Athreya and Chen ([3j), together with a number of 
other interesting theorems. We include a full proof of the new boundary Harnack principle 
because it is different from that in [3j], and ours is amenable to generalizations that will be 
the subject of a forthcoming article. 

Both techniques mentioned above — the Brownian excursion process and the boundary 
Harnack principle — are limited to Lipschitz domains and, moreover, the Lipschitz constant 
has to satisfy a certain inequality. A natural question arises whether such special Lipschitz 
domains are the largest natural family of sets where our results hold. It turns out that they 
are not. In the last section of the paper we will show that, for the two particle process. 
Too = oo, a.s., in all polyhedral domains, with arbitrary angles between the faces of the 
boundary. Unfortunately, our method cannot be easily adapted to the multiparticle case, so 
we leave this generalization as an open problem. 

For some related results on Fleming- Viot type models in smooth domains, see 16|| and 
references therein. The discrete version of the model is studied in {2]; see also references in 
that paper. 

We are grateful to Zhenqing Chen, Davar Khoshnevisan and Yuval Peres for very helpful 
advice. 
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2. Preliminaries 

For y = {y^, . . . , y^) G M'^, let \y\ denote the Euclidean norm of y and let y = {y^, . . . , y'^^^). 
We will denote the open ball with center x and radius r by B{x, r). The closure of a set A 
will be denoted A and its interior will be denoted IntA. All constants, typically denoted by 
c with or without subscript, are assumed to be strictly positive and finite. 

A function F : M'^^^ — ^ M is called Lipschitz if there exists a constant L such that 



\F{x) - F{y)\ < L\x - y\, x,y E 



Any constant L satisfying the above condition will be called a Lipschitz constant of F. 

Consider a bounded connected open set D C M*^, d > 2. We will call D a Lipschitz 
domain with Lipschitz constant L if dD can be covered by a finite number of open balls 
Bi, . . . , Bn such that for every i = 1, . . . ,n, there exists a Lipschitz function Fi : M'^^^ — > M 
with Lipschitz constant L, and an orthonormal coordinate system CSi such that 

DnB, = {{y\ ...,/) in CS, : y'^ > F,{y)} n B,. 

The following Harnack principles can be found in 

Theorem 2.1 (Harnack inequality), (a) Suppose < r < R. There exists c = c{r,R,d) 
such that if u is nonnegative and harmonic in B{0,R) C M.'^ and x,y E 5(0, r), then 

u{x) < cu{y). 

(h) Suppose that D dW'- is a domain and x,y E D can be connected by a curve 7 C -D such 
that inf^e^ dist(2;, dD) > R. There exists c = 0(7, R, d) such that if u is nonnegative and 
harmonic in D, then 

u{x) < cu{y). 

Theorem 2.2 (Boundary Harnack principle). Suppose D is a connected Lipschitz domain. 
Suppose V is open, M is compact and M G V. Then there exists a constant c = c(M, V, D) 
such that if u and v are two positive and harmonic functions on D that both vanish contin- 
uously onV r\ dD, then 

u(x) u(y) , ^ ^ 

< x,y E MnD. 



v{x) v{yy 

The next theorem is a simplified version of Theorem 1 of [l] 



Theorem 2.3. Assume that D is a Lipschitz domain. Then there exist constants tq = 
tq^D) > 0, c = c{D) < 00 and a = a{D) > 1 such that if z E dD and < r < ro then for all 
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functions u and v that are hounded, positive and harmonic on D Ci B{z,ar), and vanishing 
continuously on dD fl B{z, ar), we have 

u(x) u(y) ^ 
^<c^, x,yeDnBiz,r). 
v{x) v[y) 



Remark 2.4. Theorem 12.31 can be used to estimate the constant c{M,V, D) in Theorem 12.21 
as follows. Suppose that tq and a are as in Theorem 12.31 and we can find balls Bi{xi,ri), 
i = 1, . . . , n, and B'j{yj, p), j = 1, . . . , m, p > 0, < Tq, G dD, yj e D,M c[j- Bi{xi, ri)U 
[jj Bj{yj, p), and [j-Bi{xi,ari) C V and [jj Bj{yj,2p) C D. A simple chaining argument 
based on Theorems 12.11 and 12.31 then shows that the constant c(M, V, D) in Theorem 12.21 
depends only on n, m and D. 

I — I 

Next we recall some notation and results from ^Jj]. Fix d >2 and p > 0. Let 

(2.1) hie) = hp,d{0) = F{-p,p + d- 2; (d- I) /2;{l-cos9)/2), 
where 

F(a,6;c;x) = f;^^M^x^ N < 1, 

fc=0 

denotes the hypergeometric function and (a)^ = a{a + 1) . . . [a + k — 1), (a)o = 1. The 
function h has at least one zero in (0,7r); let 6p^d denote the smallest one. The quantity Op^d 
is strictly decreasing in p for any fixed d > 2, and strictly increasing to 7r/2 in c? for any fixed 
p > 1. In particular, if p = 2, then 

= ^^'^'^os cot^2,d = ^^^i - Therefore 6^2,2 = vr/4 and p < 2 is equivalent to 

cot 9p^d < 

For d > 2 and p > we let 6* be the angle between y and (0, . . . , 0, 1), 

Kp^d= {y eM'' ■.yj^0,0<9 < 9p^d} , 

and let O denote the axis of Kp^d- Obviously p < p' implies Kp/^d C Kp^d- We will drop the 
subscripts p and d and write K instead of Kp^d whenever there is no danger of confusion. 

The function v{x) = \x\^%{9), where h is given by (12. ip . is positive and harmonic inside 
K and continuous on K with v{x) = for x G dK. 

Let (P^,Xj) be c?- dimensional Brownian motion and for a Borel set A C M'^ define 

(2.2) Ta = mi{t>0:XteA}. 
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Lemma 2.5. Let F denote the intersection of K = Kp ^ and a hyperplane orthogonal to 
O. Let zq be the point of intersection of O with F and assume that zq G K. There exists 
c = c{p, d) such that for all zi, Z2 & O with \zq\ < \zi\ < \z2\, we have 



Proof. Let K^, be the unbounded component of K \ F and 

u{z)=¥'{Tf<T9k), zeK,. 

Then u is positive and harmonic in K^, and continuous on K^: \{F D dK), with u{z) = for 
z G dK \ F. It is easy to see that u{x) ^ as |x| — > oo. 

If I{x) = then the function u{x) = \x\'^^'^u {I{x)) is positive and harmonic in 

K = I{K^) (see Lemma 1.18 of j^). The function u vanishes continuously on dK \ I{F). 
Let K' = {1/2)K. Recall that v{x) = \x\Ph{6) is positive and harmonic inside K and 
continuous on K with v{x) = for x G dK. By the boundary Harnack principle, 

(24) u{zl v{zl 

^ ' uiz') - viz'Y 

for z, z' G K', where c depends on K and K' and does not depend on z and z'. Note that 
u{x) = Ixp^'^M {I{x)). Hence, for Zi^ Z2 E O (1 1{K'), 



\Z2 


2-d 


Z2 


-Ph{0) 




2-d 


Zl 


-Ph{0) 



u{zi) \zi\'^-'^u{I{zi)) - \zi\^-<^v{I{zi)) |zi|2-'^|zi|-P/i(0) \\zi 



The inequahty holds for all zi, Z2 E O H I{K^) (possibly with a different value of c) because 
the function u is bounded below and above on O \ I{K') by strictly positive and finite 
constants. This completes the proof of (12.30 . □ 

We will use the following estimate in the proof of Lemma 14.11 

Lemma 2.6. There exists a cone K' C K = Kp^d and a constant c = c{K, K') such that for 
X e K' and t > |a;p, 

(2-5) '"iw) ^ -^^^^'^""'^-'{W 



Proof. See or Ul. □ 
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3. A BOUNDARY HARNACK PRINCIPLE 

Let D C M.'^, d > 2, he a bounded Lipschitz domain and let A G D he a. compact set with 
IntA ^ 0. For a; G -D, define 

/(a;) = P^(T^<Taz5), 
g{x) = WTsD. 

Theorem 3.1. Assume that the Lipschitz constant L of D satisfies L < . Then there 
exists a constant c = c{A, D) such that for all x E D, 

3.1 - < ^ < c. 

c g[x) 

Remark 3.2. The condition L < ^^-j— j- is sharp. See Example 13.31 below. 

Proof of RHS of (13.11) . Since A is compact, inf^-g^ dist(x, D'^) = ci > 0. Therefore, 

inf WTqd > inf E^TaB(x,ci) = > 0. 
By the strong Markov property applied at Ta, we have for x G -D, 

E^Ta^ > c^F^Ta < Tod), 
which implies the RHS of (13.11) . □ 

Proof of LHS of (13.11) . Since Z) is a bounded Lipschitz domain with Lipschitz constant L < 
-^7=, it is easy to see that there exist p G (0, 2) and p > with the following properties. 

(i) dist(A,(9D) > 2p. 

(ii) Consider any x E D with dist(x, 91)) < p2^^. Then there exists Xq G dD and an 
orthonormal coordinate system CS = CSx^ with the following properties. The origin of CS" 
is Xo, Kp d n -B(xo, 2p) C -D n B{xq^ 2p), and x E O (that is, x belongs to the axis of Kp^^)- 
For r > and integer fc, let 

E; = {y(^W' in C S:\y- Xq\ < rtan(^p,d), |/ - x^^l < r] , 
Ek = E^-k- 

We can choose Xq and so that for some Lipschitz function F = F.^^^ : W^-^ ^ R with 
Lipschitz constant L, and all k such that 2^^^ < p, 

Z}nEfc = {(y\...,/) in C5:/>F(y)}n^fc. 

We fix X G -D with dist(x, < p2^^ and the corresponding coordinate system CS for 
the rest of the proof. 



NON-EXTINCTION OF A FLEMING- VIOT PARTICLE MODEL 



7 



Let Ek = Ek\ Ek+i and Ck = lnt{D n Ek) for k = Nq, . . . , Ni, where 

No = mm{k : 2~'' < p}, Ni = max {A; : |a;| = x'^ < 2"^"^} . 

Also let Cno-1 = Int (^D \ E^^^ and Cn,+i = lnt{D n E^^+i). 
Note that n Q = if i 7^ j, and D = Cno-i U . . . U Cn^+i- 

Let G{x,y) denote the Green function for Brownian motion killed on exiting D. Then 



(3.2) g{x)=WTao= I G{x,y)dy= ^ j 



For k = No, . . . , Ni denote by yk the midpoint of the line segment being the intersection 
of Ck with a;°'-axis in CS. In other words, {yk} = ^^-^(3/4)2-'= ^- -^^^ ^ ^'^'^ i ^^^^ ^^^^ 
j > I, k > Nq, j + k < Ni and consider the points yk and ?/a:+j- 

Let 

-^fc = c**; n Ck+i n 

and 

By Lemma [2.51 

-k \ 2-p-d 



n{yt) > ciu{yt+,) ( iptrj I = Cit<('i;t+_j)2^f^-''-'>, 



where Ci = Ci(p, rf). By scaling properties of Brownian motion, u{yk+j) = C2 = C2(p, d), that 
is, u{yk+j) depends only on p and d. We obtain 

(3.3) r (T^,,, < Tax„J > C32-^(^+'^-2), 

where C3 = C3(p, d). 
Let 

Note that v{yk+j) < 1 and v{yk) > u{yk) > C32-^^p+'^-^\ by fl33D . 

We will apply Theorem 12.21 with M = dE*^^^-^^_f,_j and \^ = -Efc+j- It follows from Remark 
12.41 that the constant C5 = c(M, V, D) may be chosen independent of k and j. The boundary 
Harnack principle implies that 

/O.N G{x,z) v{z) 
(3-4) 777 r > C5- 



G{x,yk+j) v{yk+j)' 
for 2; G -D n M. The harmonic functions G{x, ■ ) and v have zero boundary values on 
dD \ i?^3/4-)2-'=-j 5 so the inequahty (13. 4p extends to all 2; G \ E'^^^^-^^_k_j , in particular, it 
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applies to z = y^. Hence, 

(3.5) G{x,y,) ^ ^^Jp^ > c5C32-^(^+'^-^) = ce2-^^^^'-'\ 
G{x,yk+j) v{yk+j) 

Now consider the function 

hm{z) = (r^^^^ < Ten) ■ 

By the scahng properties of Brownian motion, hmiym) > C7 > for all m = Nq, . . . ,Ni. 
By the boundary Harnack principle (Theorem 12.21) applied to u{z) = G{x, z), v{z) = hm{z), 
M = Cm and V = lnt{Em~i \ -^(3/4)2-^-1)5 we have 

Gix,y) ^ G{x,ym) 

for y G C^, where Cg depends only on by Remark [2. 4[ Therefore, for y G C^, 

hmiy) 1 

(3.6) G{x, y) < cs G{x, ym) 7-7 — r < cg — G{x, y^) = cg G{x, ym) ■ 
This implies 

(3.7) f Gix, y) dy < c^G{x, y,,+,) vol(Cfc+,) < Cio 2-''^^+^^G{x, t/fc+,), 

where Cio depends only on D. 

On the other hand, by the usual Harnack inequality, 

G{x,y) > cnG{x,yk) 

for y e Bk = -B(yfe, 2"^"^), because B{yk,2^''^^) C D \ {x}. This implies that 

(3.8) / G(x, y) dy > cnG(x, y^) vol(5fc) = cu2~'"'Gix, y^), 

JCk 

where C12 does not depend on k. 

Combining (13. 5p . (13.71) and (13.81) we have 

G{x,y)dy<c,:,2^^^-^~^ j G{x,y)dy, 

where C13 = 013(1^). Fix g < 1. Since p G (0,2), we may choose j so large that ci32-'(^~^) < 
g < 1. Let Ofc = /^^ G{x^ y) dy, then 

ak+j < qak, k = No, . . . , Ni - j. 

Let N2 = min(A^i, A'^q -|- j — 1). The last inequality implies that 

A^i N2 00 N2 00 N2 

(3.9) at = afc+mjl{A:+mj<Afi} < O'fcQ'™ = C14 Ofc. 
A;=Aro fc=Afo m=0 k=N(-, m=0 fc=7Vo 
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Recall that G{x, ■ ) has zero boundary values on dD, so it is bounded by sup^g,^^^^ G{x, z) 
on the set D \ E^o- This and (13. 6p imply that sup^^^^^^ G{x,z) < ci^G{x,yNo)- We use 
(1^ to see that 



(3.10) aNo-i= G{x,y)dy < Ci5G{x,yNo)'^o\{CNo-i) 

< ci5 vo\{Gno-i)cu2^'"^ / G{x, y) dy = cieajvo- 

JCno 



Recall the definition of Nq to see that Cie depends only on D. 

The following calculation is presented in the case d > 3 only. The case d = 2 requires 
minor modifications and is left to the reader. 

Let G{x, y) denote the Green function for Brownian motion in M*^, and let G{x, y) be 
the Green function for Brownian motion in B{x, 2^^^^^). It is well known that for (i > 3, 
G{x,y) = cnlx — where cn depends on d, and G{x,y) = G{x,y) — G{x,z), for 

y e B{x, 2-^1-^) and z G dB{x, 2~^^-^). It follows that for \y - x\ < 2"^^-^, 

(3.11) Gix,y)>c^sGix,y). 

We have G{x,y) < G{x,y) for y E D, and /g^^. G{x,y)dy = cigr^. Therefore, 
(3.12) 

aNi+i= / G{x,y)dy < / G{x,y)dy < / _ G{x,y)dy = C2q2~'^^K 

Since B{x,2-^^-^) C D, 

(3.13) G{x,y)>G{x,y). 
Put yjvi+i = {x,Xd + 2~^i-5). Then by (l3TTi) and (l37[3|l . 

G{x,yN,+i) > c,sG{x,y) = C2i(2-^^)2-'^. 
Moreover, by the usual Harnack inequality, 

G{x,y) > C22G{x,yN^+i), 
for y G B{yN-^,2~''^^^^). Therefore, 



oa^i = / G{x,y)dy> / G{x,y)dy 

(3.14) -^f^^i JB(yjVi ,2-^1-2) 

> C22G(x, y^,+i) vol(5(y^,, 2-^-2)) > 023(2"^^)^-'^ ■ 2-''^' = C2,2- 
Combining fl3.12p and (13.141) . we obtain 

(3.15) oati+i < C25aNi- 
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Let C=K = C^g^i U . . . U CiVj and note that A C C^,. Let ac, = Jq^° ^{Xs&c*} ds. Then 
( K9\f . fl3J0|) and fl37[5|) imply that 



(3.16) E'Tez) < C26EVa. 

Since D is bounded, sup^g^, E'^Tg/? = C27 < 00. By the strong Markov property applied at 
the hitting time of C^:, for z & D, 

EVa < C27P"(rc. < Tqd). 

This and (13.161) yield 

(3.17) E^TeD<C28F%Tc,<TeD). 
Consider functions 

ei(^) = W{Ta<T9d), 
Uz) = r{Tc,<T9D). 

Both functions are positive and harmonic in D \ C,,,, and continuous on D\C^: with u{z) = 
v{z) = for z G dD \ C*. We apply the boundary Harnack principle with V = D \ and 
M = En^+i to see that 

/0 10N ^ ^i(yjV2 + i) 

6 (a;) 6(z/JV2+i) 

We use Remark 12.41 to see that C29 may be chosen so that it depends only on D. It follows 
from the definitions of A^o, N2 and j that for some constant C30, we have dist(?/Ar2+i, dD) > C30. 
This implies that C,i{yN2+i) = IP^'^2+^(T4 < Tqd) > C31, for some C31 depending only on D. 
We obtain from fl3.18p that ^i{x)/C,2{x) > C29C31, and this combined with (13.171) gives 

E"Taz5 < (c28/c29C3i)P^(T^ < Ta^)- 

We have proved the LHS of (13. ip for x satisfying dist(x, dD) < p2~^. 

It is easy to check that inf{/(x) : dist{x,dD) > p2~^} > and sup{(y'(x) : x G D} < 00, 
so the LHS of holds for all a; G £>. □ 

Example 3.3. The condition L < in Theorem 13.11 is sharp. To see this, note that 

for any L > ^^-j— j- there is a p > 2, such that the cone K = Kp^d is a Lipschitz domain 
with the Lipschitz constant L. Let r > be such that for every x G O, B{x^r\x\) C K. 
Then g{x) = WTqk > E''TaBix,r\x\) > Cir^jxp. Recall that f{x) = P^'(^a < Tqk) and let 
u{x) = \x\^hp^d{(^). By the boundary Harnack principle applied to / and u in a neighborhood 
of 0, f{x) < C2\x\P for X G O, \x\ < 1. Since p > 2, we cannot have f{x) > c^g^x) in a 
neighborhood of 0. The domain K is unbounded but it is easy to extend the argument to 

5(0,1). 



non-extinction of a fleming- viot particle model 11 

4. Construction of an auxiliary process from Brownian excursions 

Let f2 denote the family of all functions u : [0, oo) ^ M'^ U {5} continuous up to their 
lifetime R{(jj) = inf {t > : uj{t) = 6} and constantly equal to 6 for t > R, where 6 denotes 
the coffin state outside M'^. Let X be the canonical process on Q, i.e., Xt{u!) = uj{t) and let 
denote the distribution of Brownian motion starting from x G M.'^. As in (12.21) . for a Borel 
set A C M'^ let T4 = ini {t > : Xt e A}. Let K = Kp^a for some p > 0, and let X' denote 
the process 

Xt, for t < Tqk, 
6, otherwise, 

i.e., X' is the process X killed on exiting K. If X has the distribution P^, then X' is called 
Brownian motion in K and its distribution is denoted by P^. 

Let U denote the family of all functions u : [0, 00) ^ K U {6} such that u{0) = 0, contin- 
uous up to their lifetime R. Let denote a standard excursion law of Brownian motion in 
Kp^d starting from 0. Namely, is a nonnegative and cr-finite measure on Q such that X is 
strong Markov under with the P^^ transition probabilities and (limt^o 7^ 0) = 0. 
We have H^{n \U) =0. The existence of follows from results of 120|] and Q. 

Lemma 4.1. There exists c G (0, 00) such that 

(4.1) H^{R>t)=ct-i, t>0. 

Proof. Let =J3^ ...,0,5) G M'^ and let GK{x,y) denote the Green function for K. By 
Theorem 4.1 of [Sj, 

(4.2) H'iR>t) = c^\im^l^^j^. 

By Theorem 2.2 of Q|, which is an improvement of the boundary Harnack principle, there 
exists c{K, e) such that for all functions hi and /12 which are positive and harmonic in K 
and vanish continuously on dK., we have 

c(/f,.)-'^<^<e(A»|lM, 

h2[y) h2{x) h2{y) 

for all x,y E K n B{0, e), and limg^o c{K, e) = 1. Therefore, the limit 

h{z) 
hm — — - 

exists and belongs to (0, 00) for all functions /ii,/i2 satisfying the above assumptions. We 
apply this claim to hxiz) = Gxiz^yi) and h2{z) = \z\Ph{9), to conclude that 

lim = c G 0, cxD , 
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and 

(4.3) ffO(^>i) = „toE!l(Zk^. 



By Lemma [2.61 



C t 2 < ^ i < ct 2 . 



for t > e'^ which imphes c~^t'^ < H^{R > t) < ct'^, for t > 0. Therefore H^{R > 1) is a 
positive and finite number. 

Now mimicking the proof of Proposition 5.1 of [s], using fl4.3l) instead of fl4.2l) . we easily 
see that if {X{t), t > 0} has the distribution then for every a > the scaled process 
{y/aX{t/a), t > 0} has the distribution dP/'^H^. In particular, for every a > 

H^[R >t) = aP/'^H%R >at), t > 0, 

and putting a = 1/t we obtain (Hi!) with c = H^{R> 1). □ 

Let A denote the Lebesgue measure on = [0, oo) and let V he a Poisson point process 
on M+ X U with characteristic measure A x H^, i.e., P is a random subset of M+ x U such that 
for every pair Ai, A2 of disjoint nonrandom subsets of R+ x U, card('Pnyli) and card(Pnyl2) 
are independent random variables with Poisson distributions with means (A x H^)[Ai) and 
(A X H'^)(A2), respectively ([iS]). With probability 1, there are no two points with the same 
first coordinate, and therefore the elements of V may be unambiguously denoted by (t,et). 
Let 

Rt = inf {s > : et{s) = 6} . 
By abuse of notation, for a generic element e of U we will write 

R{e) = inf {s > : e{s) = 6}. 

Lemma 4.2. If p E (0, 2), then for every s > 0, 

-Rt < 00, a.s. 

t<s 

Proof. We use Theorem 4.6 of [18|]: if f '■ 1R+ x U ^ R_|_ is a measurable function, then 

ip{t, Ct) < 00, a.s. 

t 

iff 

/ / {fit, e) A l)dtH^{de) < 00. 

J Jr+ X U 

In particular, if ip{t,e) = i?(e)l[o,s](t), then 

Rt < 00, a.s. 

t<s 



NON-EXTINCTION OF A FLEMING- VIOT PARTICLE MODEL 13 

iff 



[ [ (i?(e) A l)dtH\de) < oo. 

J J[0,s]xU 

If we let U- = {eeU : R{e) < 1} and [/+ = {e G ?7 : R{e) > 1} then 

[[ {R{e)Al)dtH°{de) = s [ R{e)H^{de) + sH^{U+). 

JJ\0.s]xU Ju' 



'[0,s]xU 

By Lemma [4.11 



/oo noo 
H%R edt) = c j t~P/^-^dt < oo, 

because p > 0, and 

[ RdH^ = [ tH°{R edt) = c [ t- r'P/^-^dt < oo, 
Ju- Jo Jo 

because p < 2. □ 

By Lemma H72l the following process Z with values in Kp ^ U {0}, where p G (0, 2), is well 
defined. For every t > the formula 



t 



^ + Y.u<r if there exists r <t such that Y.u<r Ru < t < J2u<r ^u, 
J2u<r ^u, otherwise, 

defines a unique pair (r, s) with r > and s G [0, Rr) (in the first case). Then we define 



(4.4) Z, 



er(t), if there exists r <t such that J2u<r Ru < t < J2u<r ^u-, 
0, otherwise. 



The process Z takes values 'm. K U {0}. Let at = J2s<t t > 0. 

Lemma 4.3. The process a is a stable suhordinator with index p/2. 

Proof. The process a is increasing and has values in [0,cxd). Its paths are right-continuous 
with left limits. Note that {{t, R{et))}e<^v is a Poisson point process on R_|_ x R_|_ with 
characteristic measure A x 11, where 11 is given by 

I{{dx) = H^{R G dx) = cx-P^^-^dx, 

where the last formula follows from Lemma 14. 1[ This implies that a is a process with 
independent and stationary increments, so cr is a Levy process. Moreover cr is a subordinator. 
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since it has values in [0, oo) only. We use calculations that can be found in Section 0.5 and 
on page 73 of [7] to see that the Laplace transform of a is 

Eexp(-Aaf) = exp l-t I (1 - e^^'')Ii{dx) 







exp 



oo 



ct (1 - e-^'')x-'P'^-Hx 







= exp[ 

Therefore a is stable with index p/2. □ 

5. Construction of a Fleming- Viot process 
We recall the following informal description of a Fleming- Viot-type particle system from 



id]. Consider an open set D C and an integer N > 2. Let = {X^, . . . , X^) be 
a process with values in defined as follows. Let Xq = {x^, . . . ,x^) G -D^. Then the 
processes X^, . . . , Xf^ evolve as independent Brownian motions until the time ti when one of 
them, say, X^ hits the boundary of D. At this time one of the remaining particles is chosen 
uniformly, say, X'^, and the process X^ jumps at time ri to X!^_^. The processes X/, . . . , X^ 
continue evolving as independent Brownian motions after time Ti until the first time T2 > Ti 
when one of them hits the boundary of D. Again at the time T2 the particle which approaches 
the boundary jumps to the current location of a particle chosen uniformly at random from 
amongst the ones strictly inside D. The subsequent evolution of X proceeds in the same 
way. The above recipe defines the process X^ only for t < t^o = lim,fc^oo ^fc- There is no 
natural way to define the process Xj for t > Too- Hence, it is a natural problem to determine 
whether Too = 00, a.s. 

Theorem 5.1. There exists a constant c = c{N,d) such that if D G M."^ is a bounded 
Lipschitz domain with the Lipschitz constant L < c{N,d), then t^o = 00, a.s. Moreover, 
c{N, d) increases in N, decreases in d and 

(5.1) lim c{N,d) = c{d) ^ 



Proof. First note that Too is finite if and only if all the processes Xj;, . . . , X^^ hit dD at the 
same time, so we need to prove that this is impossible. The idea of the proof is to construct 
processes Yf^, . . . , which are easy to analyze, with values in D and such that for every 
1 < J < iV, 

(5.2) [t : Xi edD] (l{t: y/ e dD] = Aj. 

Then we will prove that that fl . . . fl = 0, a.s. 
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Recall the definition of 9p^d and Kp^^- Let -p! — 2 — 2/N and 



c{N, d) = cot 9p'^d- 



Define 



Dr^{xeD : dist(x, dD) > r} . 

Since D is Lipschitz, there exists a small r > for which the following is true. For every 
X G D\Dr there exist an orthonormal coordinate system CSx, Vx G dD and a Lipschitz 
function : M*^"^ — >■ M such that |/a; is the origin of CSx and 



Moreover, since L < c{N, d), we can choose and CS^ so that we can find a cone with 
vertex and axis passing through x which can be described in CSx as Kp^d with p < p', 
and such that 



Next we will present a very special construction of the process X^, based on a family of 
independent Brownian motions. We need this construction to show independence of processes 
Y^, Y^, . . . , Y'^ , to be constructed in a subsequent step. 

We fix j G {!,..., A^} and let denote the time of the A;-th jump of X^. We represent 
the evolution of on the interval [0, rl) as follows. We start with a family of independent 
Brownian motions W^, W^, W"^, ... in R"' starting from 0. Suppose that Xq = Xq E D\Dr/2- 
The argument needs only minor modifications if Xq G Dr/2- Let = + Xo and consider 
the cone defined as above. Let 



Inductively, for n > 1, given Xn = ^ we define the cone K^^, let PF" = W"' + Xn, and 



D n B{yx, r) C {y in CSx ■ > Fx{y)] n r). 



(5.3) 



Kp:^cL n B{yx, r) C Kp^a n -B(y^, r 



) CL'nS(y„r). 



inf{i>0: iy°G9X,J, 
a'l Ainf {t > : ^ ^r} , 




n+l 



.11 

n+l 



(j^+i Ainf > : Wl' ^ D,} 




n+l' 



n+l- 
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Then we define Tq = 0, 

n 

Tn = TlA'^ak, n>l, 

k=l 

and 

Xi = W,\ for Tfc.i < t < Tfe, A; = 0,l,2,... . 

Tliis procedure represents on tlie interval [0, t(). Strictly speaking, tI is defined in terms 
of W°,W\W\.... 

The complete construction of Xt on tlie interval [0, Too) requires that we start with a 
family {M^"'''^'"}i<j<Af,fc>o,n>o of independent Brownian motions. For fixed j and k, the sub- 
family {W^''^'"'}n>o is used to construct on the interval [r^, t^^-^), according to the recipe 
described above. The whole procedure is straightforward and elementary but tedious to 
describe so we leave the details to the reader. 

Let {Z^'^'"^''^}i<j<N,k>o,m>o,n>o be a family of independent copies of the process Z defined 
in (j4.4p . independent of {W^-''*^'"}i<j<Ar,fe>o,n>o- We will present a construction of on the 
interval [Tq,Ti) C [0,tI), based on X^ and {Z^'^'^'"'}n>o- For any fixed j,k and m, we 
can construct on [Tm,Tm+i) C [Tfc,T^+i) using X^ and {Z-^'*^'"^'"'}n>o in an analogous way. 
(Strictly speaking, we should not speak about [T^, T^+i) C [r^, t^_|_i) but about a subinterval 
of [tI, tI_^_^) constructed in a way analogous to [Tm, Tm+i)-) 

Let Tn be the isometry that maps Kp ^ onto K^^ and let = T(Z-'''^'°'"). We introduce 
a moving cone C^^ with the vertex Xj and the axis parallel to the axis of Kx„, but directed 
downwards, i.e., we put C^{x) = x — K^^ and C" = C"(X/), Tn <t < Tn+i- Let 

Si = a'[ A inf {t > : Z\ G SC^} , 
f/'^ =Xi-X(Si) + Zi(Si), t>Si, 
i?i = o'[ A inf |t > 5i : f/^^ G dK^^ } . 
By definition, S\ <T\. Then we define 




for t G [0,Si), 
for t G 



For n > 2, we define 

Sn = a'l A inf {t > Rn^i : Z'' [t - G 96^^} , 

_ -^^ _ X[Sn) + 2'" — Rn-l) , ^ > Sn, 

Rn = o'[ A inf |t > 5„ : F/") G } , 
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and 



Z^it-Rn^i), fort G [i?„_i,S'„), 

hlte[Sn,Rn). 

We continue this process until the time a" and then we put Y^ = Y^„ for t G [a", cxi). 

By construction, processes X^, . . . , and Y^, . . . , Y'^ satisfy (15.21) . Moreover, indepen- 
dence of all processes in the family {^^'^•™'"}i<j<Ar,A:>o,m>o,n>oU{Vr^'^'"}i<j<Ar,fc>o,„>o implies 
that processes Y^, . . . , Y'^ are independent. It remains to prove that, a.s., 

Ai n . . . n Ajv = 0. 

Recall that Aj = {t > : Y^ G dD}. The construction of Y^ from independent pieces of 
processes analogous to Z suggests that Aj is the range of a stable subordinator, because of 
Lemma [4.31 The matter is slightly complicated by the fact that Y^ = Y^„ for t G [(T",cri), 
and similarly for other analogous intervals. To deal with this problem, we introduce the 
following sequence of stopping times, Uq = 0, 

U; = mi{t > Uk : maj^dist(X/, dD) > r}, k> 0, 

Uk = inf{t > f/fc_i : ma5^dist(X/, dD) < r/2}, k > I. 

On each interval {Uk, Ul], sets Ai, . . . , An are independent and each one has the same dis- 
tribution as the range of a stable subordinator with index p/2. It will suffice to prove that 
for each fixed k, a.s., 



(5.4) {Uk,u*k]nA,n...nA 



N 



We use the following result of Hawkes [l7|]: The ranges of two independent stable subordi- 
nators with indices a and (3 intersect if and only if a + /3 > 1, in which case the intersection 
is stochastically equivalent to the range of a stable subordinator of index a + (3 — 1. There- 
fore, by induction, (15. 4p holds if and only if ^ — + 1 < 0. This condition holds since 
p < p' = 2 - |. □ 

Remark 5.2. Let D G M.'^ he a bounded Lipschitz domain with the Lipschitz constant L < 
Then by (15.11) we see that there exists A^o so large that L < c{N, d) for all N > Nq. 
In consequence, the Fleming- Viot-type particle process Xf in D is well defined for alH > 
provided it consists of particles with N > Nq. 

1 I — I 
Example 5.3. The proof of Theorem 1.1 in pj]] contains an error. Formula (2.1) in [IQ\ 

does not follow "by induction" from the previous statement. We will show that the error 

is irreparable in the following sense. The proof of Theorem 1.1 in [lo| is based only on two 

properties of Brownian motion — the strong Markov property and the fact the the hitting 
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time distribution of a compact set has no atoms (assuming that the starting point hes outside 
the set). Hence, if some version of that argument were true, it would apply to almost all non- 
trivial examples of Markov processes with continuous time, and in particular to all diffusions. 



However we may find a diffusion for which the analogue of Theorem 1.1 in 10| is false. Let 
Xf be the diffusion on [0, oo), started at Xq = 1 and satisfying the SDE 

dXt = dWt - -^dt. 

2Xt 

We make absorbing so that it can play the role of the boundary for the domain D = (0, oo). 
Notice that although Xf is not a Bessel process, as we have reversed the drift term, it scales 
in the same way. That is, for a > 0, aXfa-^ is a diffusion satisfying the same SDE, but 
started at a. Let = [Y^'^ , 1^/'^) , i = 1 . . . oo, be a double sequence of independent copies 
of Xt, and set 

ai = inf {t > : F/'^ A F/'^ = O} , 

Now, construct a two-particle Fleming- Viot type process = {Xl,X^) as follows. First 
let Ti = (Ji and set = Y^^ for t E [0, ti). At ti one of the particles hits the boundary 
and jumps to = ai. To continue the process we use the scaling property of Y^ and set 
Xi = ^lYf 2 for t E [ti,T2) where T2 = Ti + ^i^o-2- At T2 a second particle hits the 
boundary and jumps, this time to ^2 = 02^1, and we continue the process in the same way 
by setting 

i i 

i=i i=i 
Xt = ^iYJ^_^^^^-2, for t E [n, Ti+i) . 

Then Xj evolves as two independent copies of Xt with Fleming- Viot type jumps when a 
particle hits the boundary. The process X^ is well defined up until t^o and if the analogue 
of 10 



Theorem 1.1] were to hold for this process we would have Too = 00 almost surely. 
In fact the opposite is true. We will show now that Etoo < 00 and hence Too < 00 almost 
surely. To do this it will be sufficient to show E (ai^) < 1 and Eai < 00. Let f {x,y) = 
x'^ + y"^ — x^y"^ and notice f{x,x) = f{x,0) = f{0,x) = x^. We may check using Ito's 
formula that / {Yt/^^.,Yt/^^.) is a positive local martingale and hence a supermartingale. By 
the optional stopping theorem 



E {a,') = Ef [Y^YXf) < Ef W'\ro^'^) = 1. 
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Furthermore, ai is not almost surely constant and so by Jensen's inequality 



E (ai^) < v^E (a/) = 1. 

We may use Ito's formula again to show that Xt^ + 4t is a local martingale and so by the 
optional stopping theorem again we have that E (cti) < j. 

By independence of the Y* processes we have that E (^j^) = E (ai^)* and so 



oo ^ oo 

Etoo = ^E < T [ai^Y < oo. 



4 

j=l j=0 

6. Hitting probabilities of compact sets 

This section is devoted to a technical estimate needed in the proof of Theorem 17.11 Recall 
definitions of Dr and = {X^, . . . , X^^). 

Lemma 6.1. Fix N > 2 and let D dW^ he a hounded Lipschitz domain with the Lipschitz 
constant L < c{N, d) . 

(i) For any fixed k G {1, • • • ,N}, and for every r > such that IntZ^^ 7^ 0; there exist 
c > and t > such that for all x G , 

(Xf G D,) > c. 

(a) For every r > such that IntD,,. ^ 0, there exist c > and t > such that for all 
xGD^, 

(Xt G D^) > c. 

Proof, (i) Fix r > such that IntD^ 7^ 0- Recall that notation such as To^, Tgo, etc. refers 
to hitting times by Brownian motion. By Theorem 13. II there exists Cq = Co(r) such that for 
all x G -D, 

(6.1) F^{Tn,.<T9n)>CoE^Ten. 
Fix k and let T^^" = inf {t > : X/^ G Dj, and 

Y, = X'{tAT^'). 



Define Tn = and 



Let Mo = and 



and 



Tn+i = inf |t > r„ : \im G dD^ A T^" . 



Mn = — l{y(r„)6Z)r} - Tn, n> 1, 

Co 



J^n = 0"(Xt, t <Tn). 
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It is easy to see that ETn < oo so E\Mn\ < oo. For x = {xi,X2, ■ ■ ■ ,X]\f) G with Xk ^ D^, 
(M„+i - M„ I 

= — E'' (l{y(T„+i)GD,}(l{y{T„)^D,} + l{y(T„)eD,}) - l{y{T„)eA-} I ^n) - E'' (T^+i - \ Tn) 

Co 

= — E'' (l{y(T„+i)eD,}l{y{T„)^D,} + l{y(T„)GD,} - l{y(T„)eD4 | ^n) - E'' - T„ | Tn) 

Co 

= li|5.(^^)^^^}P- (r(r„+i) G z^, I - E" (T„+i - T„ I Tn) ■ 

Co 

We have on the event {Y{Tn) ^ Dr}, 

E" (M„+i - M„ I > lp^'(^") (T(AO < Tqd) - E^'(^")Tai5 > 0, 

Co 

by (16. ip . On the event {y(T„) G -Dr}, we have T„+i = T„, 1t„+i £ -Dr, and so 

E'^(M„+i - M„ I = 0. 
Combining the last two formulas, we conclude that {M„} is a submartingale with respect to 

Define 

S = inf {j : T,- > 1} A inf { j : F^^^, G } . 
Fix an X G and consider two cases. First, we may have 

P^(S = inf{j:Fr, G^r})>l/2. 

In this case, 

(6.2) F^iT^' < 1) > 1/2. 
The second case is when 

P"(5 = inf{j:FT, eA.})<l/2. 

In this case, F^{S > 1) > 1/2, so W^Ts > 1/2. The submartingale M„ is bounded above by 
1/co so we can apply the optional stopping theorem to obtain 

WMs > WMo = 0. 

Hence 

(6.3) P^ {Yt, G Dr) > CoE^Ts > Co/2. 
We will show that for some to, 

(6.4) P^ (t^; < to) > Co/4. 
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If Ts > Sq for some Sq > 1 then Xj^ must not hit U dD for t G (1, Sq). The probabihty of 
this event is bounded above by the probabihty of the event that Brownian motion starting 
from win not leave the baU B{Xi, 2 diam(D)) for sq — I units of time. The last probability 
is Ci < 1, depending on Sq > 1, but not depending on X^. By the Markov property, 

sup P'' {Ts > So) < ci < 1. 

Applying the Markov property repeatedly at times sq, 2so, . . . , we obtain for any x G , 

F^{Ts > nso) < c^ 
We choose n so large that c" < Co/4 and let to = "'^-Sq. Then for x G , 

(6.5) F^{Ts > to) < Co/4. 
We use (16. 3p and (16.51) to see that 

Co/2 < {Yt, G Dr) 

= {Yt, G Dr, Ts > to) + (Ft, G D^, Ts < to) 
<¥-iTs>to)+¥- (t^" <to) 

< Co/4 + p^ (r£ < to) . 

This implies (16.41) . We combine the two cases, that is, (16. 2p and (16. 4p . to see that for some 
ti < oo and C2, for all x G -D^, 

(6.6) p^ (r£ < ti) > C2. 

Let ri be such that < r < ri and IntD,.-^ ^ 0. Let t2 and C3 be such that (16.60 holds 
with ri,t2 and C3 in place of r,ti and C2, i.e., 

(6.7) P'' (t£; < t,) > C3. 

Let r2 = (ri — r)/2 and pi = P° (T5B(o,r2) ^ ^2) > 0. By translation invariance of Brownian 
motion, pi = P^ {TdB{y,r2) > ^2) for every y. If the process hits Dr^ before time ^2 and 
then stays in the ball B{X''(T^^^),r2) for at least ^2 units of time then X^ will be inside Dr 
at time ^2- By the strong Markov property apphed at the stopping time T^^^, we obtain, 
using ([621), for all x G D^, 

(6.8) P^(Xf^ G D,) > piP'' (r^;^ < ^2) > P1C3 > 0. 
This completes the proof of part (i) of the lemma. 
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(ii) Recall that r > is fixed and such that IntDr 7^ 0. Let and r4 be such that 
< r < ra < r4 and IntDr^ 7^ 0. Let = min(r3 — r, r4 — r^)/2. We choose ^3 and C4 so 
that (16.81) can be applied with in place of r, 

P"(X^^ e Dr,) > C4 > 0. 

Let P2 = infyeD P^l^az) < ^3) and note that p2 > 0. Let ps = {TdB(y,r^^) > Sta) > and 
note that ps does not depend on y. 

Let A be the intersection of the following events. 

(a) The process is in Dr^ at time ts, and it stays in B{X^^, r^) for all t G [ts, St^]. 

(b) For every j = 2, . . . , N, the process jumps at a time Sj G [ts, 2^3] to X^., and then 
stays in the ball B{Xi^,r^) = rs) for all t G [sj, sj + 21^]. 

By the strong Markov property and the definition of the process X, the probability of A 
is bounded below by C5 = C4P3{p2{l/{N — l))p3)^^^. If A occurs then Xa^g G . Hence, 
for every x G Z^^, 

P'^ (X3t3 G D^) > C5 > 0. 
This proves part (ii) of the lemma. □ 

7. Stationary distribution for the particle system 

The two theorems proved in this section generalize the analogous results in pii], where 
the proofs were given only for domains satisfying the internal ball condition. 

Theorem 7.1. Suppose that D G Mf^ is a bounded Lipschitz domain with the Lipschitz 
constant L < c{N,d), where c{N,d) is as in Theorem \5.1[ Then there exists a unique 
stationary probability distribution for Xj. The process X^ converges to its stationary 
distribution exponentially fast, i.e., there exists A > such that for every A C D'^ , 

(7.1) lim e^* sup |P^ (X^ e A) - M^iA) \ = 0. 

Proof. We have shown in Lemma [6.11 (ii) that for any r > 0, with probability higher than 
Po = Pair) > 0, the process X^ can reach the compact set within to > units of time. 
This and the strong Markov property applied at times 2to, 4to, 6to, . . . show that the hitting 
time of is stochastically bounded by an exponential random variable with the expectation 
independent of the starting point of X^. Since the transition densities p^(x, y) for X^ are 
bounded below by the densities for the Brownian motion killed at the exit time from , 
we see that p^(x, y) > Ci > for x, y G . Fix arbitrarily small s > and consider the 
"skeleton" {X„s}n>o- The properties listed in this paragraph imply that the skeleton has a 
stationary probability distribution and that it converges to that distribution exponentially 
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fast, i.e., (17. ip holds for the skeleton, by Theorem 2.1 in [l5^ or Theorem 16.0.2 (ii) and 



of 



VI 



21| . See the proof of Proposition 1.2 in [9j for an argument showing how to pass from 
the the statement of uniform ergodicity for the skeleton to the analogous statement for the 
continuous process t — X^. We sketch this argument here. Take any e > and find ti = nis 
such that 

(7.2) e^* sup (Xt eA)- M^iA)\ < e 

holds for t > ti of the form t = ns. Consider an arbitrary t2 > ti, not necessarily of the 
form ns. Let m be the integer part of t2/s and let u = t2 — ms. Note that m > rii. Since 
(17. 2p holds for t = ms, the semigroup property applied at time u shows that (17. 2p holds also 
at time ^2- D 

Theorem 7.2. Suppose that D is a bounded Lipschitz domain with the Lipschitz constant 
L < . For N > Nq (see Remark 15.^) let Xj^ be the stationary empirical measure. 
Let if be the first eigenfunction for Laplacian in D with the Dirichlet boundary conditions, 
normalized so that Jj^ip = 1. Then the sequence of random measures Xj^, N > Nq, converges 
as N oo to the (non-random) measure with the density ip, in the sense of weak convergence 
of random measures. 

Proof. Recall processes defined in the proof of Theorem 15.11 By construction, we have 
dist(y/, dD) < dist(X/, dD), for all j and t. 

It is elementary to see that the process Z constructed in Section H] has the property that 

1 /•* 

hmhmsup- / l{dist{Zs,dD)<r}ds = 0, a.s. 

rlO t^oo t Jq 

In view of the construction of from independent copies of Z, we also have, for every j, 



1 



limlimsup- / l|dist(y/,aD)<r}^^ = 0' a.s. 



Hence, for every j, 



1 /•* 

hmhmsup- / l{di.t(x|,9D)<r}^* = 0' a.s. 

^J-'J t^oo 1^ Jo 

This implies that for every pi > 0, one can find r > so small that if X has the stationary 
measure Ai^ then for every t, P(X/ ^ Dr) < pi- It follows that for any A^, the mean measure 
EXj^ of the compact set is not less than 1 — pi. Hence, the mean measures EXj^ are 
tight in D. Lemma 3.2.7, p. 32, of [13i] implies that the sequence of random measures Xj^ is 
tight and so it contains a convergent subsequence. 
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One can complete the proof of the claim that the random measures converge as 
^ oo to the measure with the density ip exactly as in the proof of Theorem 1.4 in 101 ]. 
starting on line 9 of page 699. □ 



8. Polyhedral domains 

In this section we show that the Lipschitz constant c (A^, d) in Theorem 15.11 is not sharp, 
that is, Too = oo, a.s., in some Lipschtz domains with arbitrarily large Lipschitz constant. 
Specifically, we will demonstrate the existence of the two particle process for all times in 
arbitrary polyhedral domains. Unfortunately, our method cannot be easily adapted to the 
multiparticle case, so we leave this generalization as an open problem. 

Definition 8.1. We say an open set D C M*^ is a polyhedral domain if there exist simplicial 
complexes K. D OK, such that D = \1C\ and dD = \d]C\. 

For the remainder of this section we will assume that D = Int |/C| is a polyhedral domain. 
Let = {Xl,Xf) be a Fleming- Viot process in D and define jump times Tj as before. We 
will show: 

Theorem 8.2. If D is a polyhedral domain and Xj = {Xl,X^) is a Fleming-Viot process 
with jump times ti then Ti ^ oo as i ^ oo almost surely. 

As Xf is a CADLAG process we have X^^ = X^^ for each i G N, so we may define a sequence 
of jump points = X^^ = X^^. Since D is compact, C,i has at least one limit point in D. To 
prove Theorem 18.21 we will examine the behavior of Xj when both particles are close to a 
limit point of C,i and, assuming that r^o < oo, arrive at a contradiction. 

First we will show that if t G [rj,rj_|_i) then X^ cannot stray too far from 

Lemma 8.3. Set = \\X^ Vt = W^t - fort E [r^, Ti+i). Ifr^ < oo then V} 

and — > as t ^ Too- 

Proof. It suffices to consider only Vf^. Notice that V^^ is a (i-dimensional Bessel process 
(Bes (d), for short), reset to at each Tj. So setting AV^^ = V^. we may extract a Brownian 
motion 

{im:Ti<t} * 

Consider e > 0. We will count the number of upcrossings of the interval within a 

short time interval \t,t + 5], where 5 = £^/(4 {d — 1)). Consider times t<s'<s<t + 6 
where V/ > e and s' = sup {s < s : ^/ = f }• Notice as only jumps downwards there is 
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no i G N such that s' < Ti < s. We have 



e d- 1 

^ e d — 1 e 

- 2 ~ 4(d-l) e ~ 4' 

So on a short time interval, each upcrossing of by corresponds to an oscillation 

of I by W . As is a Brownian motion, with probability 1, makes only finitely many 
upcrossings of in a given time interval + 5]. If t^o < oo, we may find n G N with 

Tn > Too — 5- So iiV^>e for some Tn < < t < Tj+i then as F"*^ is reset to at there 
must be an upcrossing of in the interval [rj,rj+i) C [r„,r„ + 5]. So V^"*^ > e in only 

finitely many intervals [r^, Tj+i) and, as e is arbitrary, V/ ^ as t ^ Too- □ 

Corollary 8.4. If Too < oo then the sequence C,i has no limit point G D. 

Proof. Fix X E D. As D is open, there exists some e with B {x, 2e) G D. If G 5 {x,e) 
then, as both particles follow continuous paths until one exits D, we must have V^^ \/V^^>e 
for some t G [r^, Tj+i). So if VJ^ as t ^ Too then C,i E B (x, e) for only finitely many 
i. As X is arbitrary we see that so long as VJ^, V^^ ^ as t — Tqo, C,i can have no limit point 
in D. □ 

Corollary 18.41 is similar to a result in [l^ (Step 1 of Theorem 7). In that paper, a system 
consisting of an arbitrary number of particles is considered, but the boundary dD is assumed 
to be smooth. 

It is convenient at this point to introduce some notation that will allow us to consider the 
behavior of when it is close to the boundary of a simplicial complex. Let a be a fc-simplex 
with vertices {vq, . . . ,Vk}, that is 



a 



kvi ■ Ao,...,Afc > 0, y^Aj = 1 



i=0 i=0 

Then define the interior of a 

k k 

>^iVi : Ao,...,Afc > 0, y^A. 



a 



i=0 



and the span of a to be the subspace 

r k k 

5<x = < ^ XiVi : ^ Xi 



i=0 i=0 
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For two simplices (7i , (72 G /C we write ai < 02 if 0\ is a face of and a\ < if 0\ is a 
proper face of 02- We name the star of a simplex a to be the set 

St [a) = {(Tl e /C : cTi > a} 

and define the neighborhood of a as 

jV ((j) = |x e : X e (Tl for some (Ti > (t| . 

Given simphces o < 0\ name the vertices of cr and Ui, {vq, • . . , Vjt} and {vq, ...,Vn} 
respectively. Define the wedges 

{n n \ 

AiWi : Afc+i, . . . , An > 0, Aj = 1 L 

>V((t)= U W((7,(7i). 

o-ie5t(o-) 

Notice that M {a) C W (a) and that J\f (cr) is open with respect to the subspace topology 
of D. Notice also that W (cr) is a product space 

W ((t) = C ((t) X 5,, 

where the cone C (cr) is the projection of W (cr) onto S^. 

Now, consider cr e 9/C and suppose there exists a subsequence Ci™ Coo £ cr- As ^00 £ cr c 
jV (cr) and jV (cr) is open in D we may assume without loss of generality that G N (cr) 
for each n. So consider Xj started at (Ci„)Cin) time r^^ and stopped at the first time 
T > where one of X}, Xf exits A^(cr). Of course, as Af (a) C W (cr) n !D, this has the 
same distribution as a Fleming- Viot process in W (cr) started and stopped in the same way. 

So, let and be the probabihty measure and expectation operator associated with 
a Fleming- Viot process in W (cr) started at Xq = (x, x) . The Sa and components are 
not quite independent as they have the same jumps, but P^. allows a partial factorization as 
follows. 

Lemma 8.5. //Xt is a Fleming-Viot process in W (cr) then there is a well defined decom- 
position Xt = Yt + Zt with Yt = (F/, F/) e C {af , Zt = {Z^, Zf) e with the following 
properties 

• Yf is a Fleming-Viot process in C (cr); 

• there exists a Brownian motion Zf in S„ ( not adapted to the filtration of X^ ), inde- 
pendent ofYt, such that for each i & N we have Z^-. = Q with Q — Z^. — Z^.. 

Proof. Obviously, as C (a) C 5^, the factorization X^ = Yj + Zf is unique. Further, on each 
interval [ri,ri+i), the processes F/, Y^, Z\ and Zl evolve as independent Brownian motions 
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on and Sa respectively. So as is a subspace and has no boundary, jumps when 
and only when Yf hits dC{a), and so is indeed a Fleming- Viot process on C {a). 

Now for each i e N only one of X}, Xf has a discontinuity at Tj+i, so there is a well 
defined sequence of random variables Ji E {1,2} such that X/* is continuous on the closed 
interval [Ti,Ti+i] and we may define a continuous process 

Zt = z^', t e [Ti,Ti+i] . 

Then = Q for every i and it remains to show that Zt is a Brownian motion independent 
of Y^. Of course Zf is only defined up to Too- But we may continue Zf after r^o with an 
independent Brownian motion if necessary. 

Now as Zt follows either Zj or Z^ then the quadratic variation {Z'^^ = tl and, by Levy's 
characterization, we need only check that Zt is a martingale with respect to its own natural 
filtration and is independent of Y^. Furthermore, although Zt is not adapted to X^, for each 
Tj, the path Z 1^^^ is measurable with respect to Xj^^^. Therefore, by the strong Markov 
property, it is sufficient to consider only intervals [rj,rj+i). 

In fact it suffices to consider only the first time interval [0, ri). Let Xt be a Fleming- Viot 
process started at .^o ^ W (c) and stopped at ti. Then the left limit process is a pair of 
independent Brownian motions stopped at r = t{ . Set J = Jq and we have = X^ e W (cr) 
and X^--^ G dW (a). 

So set Xf = Yf + Tit as in the statement of the lemma and let ^ and Tf be 
the natural filtrations of Y, Z and Z respectively. Set Co = Ci = to be the T^- 
measurable Z-components of and respectively. Thus, r is a stopping time of and 
J is measurable with respect to . Now crucially Y and Z are independent processes so 
for t < T we have 

Thus Z is a martingale, and hence a Brownian motion, with respect to the filtration 
Qt — V Tf. Therefore Z is independent of C and is a Brownian motion with 
respect to its own natural filtration Tf C Qf 

□ 

Now Y^ is a process in a cone and if converges to some point in a then Y^ must converge 
to the apex of C (cr). Our next step is to show that this cannot be the case. 

Lemma 8.6. // Y^ is a Fleming- Viot process in a cone C C M'' then, with probability one, 
Yt does not converge to (0,0). 
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To prove this we will need to consider the angular components, $^ = n j,, , of Y. We will 

\rt II 

recall briefly^ some facts about spherical Brownian motion. We will omit details, which can 
be found in [23, Chapter 8], particularly Example 8.5.8. 

Let Bt be a Brownian motion on R*^, let the unit sphere be denoted 

§^-1 = {x e : ||x|| = 1} , 

and define the map (f) : R'^\ {0} S^^^^ by (x) = 
Now let $t = (j){Bt). Applying Ito's formula, 

d^t = jr^ (I - $t$7) dBt - ^~ dt. 
\\Bt\\ ^ * * ^ * 2\\Btf 

Note we are interpreting $t as a column vector so ^t^J is a square matrix. Now define a 
differential operator A : R) C° (§'^-\ R) by 

= H - E ^ ) - ^ E ^- 

Applying Ito's formula again, we see that / — ^^^^dt is a local martingale for each 

We may extend this to functions of two Brownian motions by defining A^,A^ by 

A'f{x,y) = A{f{;y)){x), 

A'f{x,y) = A{f{xr)){y). 

Then by a similar application of Ito's formula, if B^ and B^ are independent Brownian 
motions and $J = (5^^), = (5^^), = {^j, ^f), then 

is a local martingale. 

Now apply a time change to as follows. If ait) = inf |s G : /q* > t}, 

then 0^ = is a Markov diffusion on E>'^~^ with generator A. Let Pg^'^^ and Eg^'^^ be the 
probability measure and expectation operator associated with two independent copies of Gf 
started at 6i and 62 G S"^"^ respectively. 

Lemma 8.7. Let U be an open subset ofS'^^^ and set 

Tf = inf {t G M : 0^^ G dU} , 
= inf {t G M : G dU} , 
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Then G (f/2,M) and A^h^ = -A^h^ > 0. 

Proof. The process (Gj, Gj) is a Markov diffusion with generator A"^ + A?, so by Dynkin's 
formula A^h^ + ^^/i^ = and it remains to show that A^h^ > 0. 
By definition of the Markov generator 

A'h'' (^1, 02) = lim ^ (e^^^ (/i^ (G^, ^2)) - {0u O2)) 

= lim i (e^^^ (p,«' '^^ [Tf < T,^] ) - P^^^ [Tf < T,^] ) . 

But Eg^'^^ (P§*'^' (■)) is the probability measure associated with the process (Gj^^,G^), 
s > 0, obtained by giving Gi a headstart. So we have 

K"'" [^^i"" < ^2""]) > K"'" [Ti - t < T^] - Pg^''^ [Tf < t] 

> Pg^'^^ [Ti^ < T2^] - Pj'^^ [Ti^ < t] 

and, since jPg^'^^ [T^ < t) — as t — > 0, we may pass to the limit, and we see that 
A^h^ {9\9^)>0. □ 

We are ready to prove Lemma 18.61 

Proof of Lemma \8. 6[ Set C = {An : A G M^, u G U} for some open subset U C S'^^^ and let 
be a Fleming- Viot process in C. 
We deal first with the special case when d = 1, in which case either C = R and there is 
nothing to prove or C = M^. If C = R"*" then is a 2-dimensional Brownian motion in the 
quarter plane with jumps (y, 0) 1^ {y,y) or (0,y) 1— » {y,y) whenever the process exits the 
first quadrant. As these jumps only increase ||Yt|| then ||Yt|| dominates a Bes (2) process 
and Y^ does not converge to 0. 



For d > 2 define a function 



log a d = 2, 

LL_ if > 3, 



and define processes 

Ht = {^\, $2) , St = Ml + {Ml - Ml) Ht. 

Now, fi is harmonic on R*^ and it will be key to our argument that Ml and Mf are both 
local martingales except when Y^ jumps. We say a Y^-adatapted process Rt is a martingale 
between jumps if Rt — J2{ieN-Ti<t} {^n — -^r.-j is ^ continuous local martingale. The process 
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St is a convex combination of and M^, so if both Y^^ and Y^^ converge to the origin, then 
St converges to — oo. Notice also that if Y^ approaches dC then Ht —>■ 1 and so St —>■ M^. 
Similarly, if Y^^ approaches the boundary then St — > M^. So St is continuous. 
Set 

Jo V \\B^\\ \\B^\\ J 

By (IS.ip Nt is a martingale between jumps. We may check that the cross variation terms 
(M\ $1)^ = (M^, $2) = and so, as Ht is a function of <I>J and we have (M\ H)^ = 
{M^, H)^ = and for s E [xj, Xj+i) we may calculate 

= Sr, + / (1 - ^t) ^^t' + / Ht dM^ + I (M2 - M/) rfi/i 
= Sr,+ (1 - Ht) dMl + Ht dMf + / {Mf - Ml) dNt 

J Ti J Ti J Ti 

Therefore 5*8 — {M^ — M}) i ^^-f^ + ) is a martingale between jumps. 

Now from Lemma [8.71 we have A^h^ = —J?hF > and so 



.m\ \\Bt.^ 

But /i is an increasing function of the norm || ■ ||, so for Tj < Si < S2 < Tj+i, 

Therefore St is a continuous local submartingale and it cannot converge to — oo. Thus Yt 
does not converge to (0,0). □ 

Corollary 8.8. //X^ is a Fleming- Viot process in a polyhedral domain D then with proba- 
bility one the sequence of jump points does not converge to any ^oo £ dD as i oo. 

Proof. First, for a G dJC, let be the event that C,i C,oo for some ^oo £ o" and assume 
without loss of generality that Oe a. Set 

i^'^ = n [xi eAf{a);t> n, j = 1,2] . 

Then, as A/" (a) is open in D, from Lemma 18. 3[ increases to F'^ up to an event of 



probability 0. By the strong Markov property and Lemma [8.61 

P {Fn = (Vt ^ (0,0) n [X{ eAf{a);t> Ti, J = 1, 2]) = 0. 
So as d)C is a finite set of simplices we have P [3^oo £ dD s.t. ^ ^00 as i — 00] = 0. □ 
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To complete the proof of Theorem 18.21 we consider the set 




a G /C : there exists a subsequence 



^ G 0" as n — > oo 



} 



It is easy to check that the event {a G L} is X-measurable. We say a is a local maximum of 
L if L n St (a) = {a}. Of course any non-empty subset of a finite lattice contains at least 
one local maximum, and L is non empty by compactness of D. We will prove Theorem 18.21 
by showing that for each cr G /C the event that Too < oo and cr is a local maximum of L has 
probability 0. 

Proof of Theorem \8.2 . Fix a G d)C, and note that TV (cr) \cr is non empty. We show first that 
if has a limit point in a and r^o < oo, then has a second limit point in A/" (a) \a. 

First suppose that a = {v} is a vertex of /C and f is a limit point of ^j. By Corollary 18. 8[ 
the sequence does not converge to f as i oo, so we may choose e > such that 
B {v,e) n D C A/'((t) and that — "^11 > ^ infinitely often. If this is the case then there 
are infinitely many pairs such that C,i„ G B{v,e) and ^ B{v,e). But from 

Lemma [8731 we have — ^ as i — > oo hence — f || ^ e as i ^ oo. Therefore, as 
dB {v, e) is compact, must have some limit point in dB {v,e) H D C M (a) \ {v}. 

If 0" is a /c-simplex for < k < d then for each x G o" there exists e > such that 
B {x,2e) n D C A/'(cr). We will consider upcrossings of the interval [e,2e] by — 
Define sequences ^„, jn G N U {oo} and T„, G M U {oo} by: jo = 0, 



Then we put = sup {n G N : jn < oo} to be the number of upcrossings. 

Note that B {x,2e) HD C A/" (a) and so X(f_|_T-,„)AT„ is a Fleming- Viot process in W (a) 
started at and stopped on exiting B{x,2e). So we may consider P|'" and factorize 

X.t = Yf + Zf as in Lemma [8751 For t G [rj,^,?7„], the process is measurable with respect 



Jn 



n 



inf {i > in ■■ B {x,e)} , 

inf {j > in ■ iB{x,2e)}, 

inf [t > Ti„ : X] ^B {x, 2e) or ^ B {x, 2e)] , 

sup {Ti : Ti < Tn} . 



to X L , which is distributed according to P^" . Hence Z 



is a Brownian motion 



in Sfj with respect to its own natural filtration. 



Recall Zt-. = (Ci, 0) and set 




32 MARIUSZ BIENIEK, KRZYSZTOF BURDZY AND SAM FINCH 

Then Vt is dominated by a Bes (d) process reset to zero at times Xj^ . So arguing as in the proof 
of Lemma [873| if < cxd and the number of upcrossings = oo, then r^^ < < oo for each 
n G N, and as ^ oo. But r]^ = sup^ {r^ : < r„}, hence X^,^ = {^k„,^kj for some 
kn and either ||x — X^^^ || = 2e or ||x — X^^ || = 2£. So if Too < oo and = oo, we must 
have — a^ll ^ 25 as n — oo and so has a hmit point ^oo ^ dB (x, 2e) H D C A/" (a). 
But Vt ^ as t Too with probabihty one, so we cannot have ^oo G o" and we must have 
^oo G AT (a) \a. 

Now let Q"^ be a countably dense subset of a and suppose has some hmit point x G a. 
By Corollary 18. 8[ C,i does not converge to x as i — > oo and we may choose some rational 
e > such that B{x,3e) H D C M (a) and — x|| > 3e infinitely often. Now choose 
q G Q'^ n B {x,e) and notice that — g|| makes infinitely many upcrossings of the interval 
[e,2£:]. If Too < oo then as Q'^ is countable, with probability one we may find some limit 
point ^oo e A/" (cr) \cr. 

Recall the definition of the set L. As L is nonempty there must exist some local maxi- 
mum a. However if Too < oo then, by Corollary 18.41 we have L C dlC. We have just shown 
that if Too < oo then L has no local maximum in d)C. Hence we must have Too = oo. □ 
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